Spontaneous decay of an excited atom placed near a rectangular plate 
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Using the Born expansion of the Green tensor, we consider the spontaneous decay rate of an excited 
atom placed in the vicinity of a rectangular plate. We discuss the limitations of the commonly used 
simplifying assumption that the plate extends to infinity in the lateral directions and examine the 
effects of the atomic dipole moment orientation, atomic position, and plate boundary and thickness 
on the atomic decay rate. In particular, it is shown that in the boundary region, the spontaneous 
decay rate can be strongly modified. 
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PACS numbers: 42.60.Da, 32.80.-t, 42.50.Nn, 42.50.Pq 



The ability to control the spontaneous decay process 
holds the key to powerful applications in micro- and 
nano-optical devices. Effective control can be achieved by 
tailoring the environment surrounding the emitters. In 
theoretical analysis of surrounding environment of differ- 
ent geometries, the most interesting ones being of the res- 
onator type, the boundary conditions are typically taken 
into account only in directions in which the electromag- 
netic field is confined or affected the most. For example, 
in a planar configuration, only the boundary conditions 
in the normal direction are taken into account while those 
in the lateral directions are neglected (see, e.g., Ref. [l[). 
In a cylindrical configuration that extends to infinity, the 
reverse is true 0. Under appropriate conditions, these 
approximations are generally valid. However, as the sizes 
of devices decrease and fall in the micro- and nano-meter 
ranges as in the current trend of miniaturization, it is 
clearly of great importance to keep track of the effects 
of all boundaries. One way to calculate the spontaneous 
decay rate in an arbitrary geometry is to directly solve 
the Maxwell equations using the finite difference time do- 
main method Q. This method, which relies entirely on 
numerical computation, is not without weaknesses. It re- 
quires that the whole computational domain be gridded, 
leading to very large computational domains in cases of 
extended geometries, or in cases where the field values at 
some distance are required. All curved surfaces must be 
modelled by a stair-step approximation, which can intro- 
duce errors in the results. Additionally, the discretization 
in time may be a source of errors in the longitudinal field 
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Here we employ an approach that, in a sense, com- 
bines analytical and numerical calculations, thereby sig- 
nificantly reducing the numerical computation workload. 
This approach relies on first writing the atomic decay 
rate in terms of the Green tensor characterizing the sur- 
rounding media [J, Q ■ Although this formula holds for 
arbitrary boundary conditions, exact analytical evalua- 
tion of the Green tensors for realistic, finite-size systems 
can be very cumbersome or even prohibitive. Following 
@ and 0] , where the atom-body van der Waals force and 
the local-field correction, respectively, have been consid- 
ered, we circumvent the task of an exact calculation of 



the Green tensor by writing it in terms of a Born se- 
ries and restrict ourselves to leading-order terms. The 
boundary conditions enter the theory only via the inte- 
gral limits. This approach is universal in the sense that 
it works for an arbitrary geometry of the surrounding 
media, and can be used to evaluate any characteristics 
of the matter-electromagnetic field interaction express- 
ible in terms of the Green tensor. In this paper, we are 
concerned mostly with the spontaneous decay rate of an 
excited atom placed near a rectangular plate. Our aim 
is twofold: first, we compare our results with those for 
an infinitely extended plate in order to establish in a 
quantitative way the conditions under which the approx- 
imation of an infinitely extended plate is valid; second, 
we examine the effects brought about by the presence of 
the boundaries in the lateral directions. 

The (classical) Green tensor of an arbitrary dispersing 
and absorbing body satisfies the equation 



HG(r,r',uj) = S{r-T')I, 

1 uj 2 

H(r) = V x — -V x — - e(r,w), 

fi[r, to) & 
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(J-unit tensor) together with the boundary condition 
at infinity, where e(r, w) |/i(r, oj)] is the frequency- 
and space-dependent complex permittivity (permeabil- 
ity) which obeys the Kramers-Kronig relations. 

Decomposing the permittivity and permeability as 

e(r,uj) = £(r,cj) + x e (r,cj), fi(r,oj) = fl(r,uj) + x fJ ,(r,uj), 

(3) 

and assuming that the solution G(r, r', iS) to the equation 

H{y)G(y, r', oj) = 5(r — r')I, where H is defined as in Eq. 
@ with e instead of e and fx instead of /i, is known, the 
Green tensor G can be written as 

G(r, r', oj) = G(r, r', oj) + G'(r, r', oj). (4) 

Substituting Eq. ^ into Eq. |T|) and using the identity 
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(M + Xn) 1 = M 1 S^CW/- 4 )*! it can be found that 
H(r)G'(r, r', w) = ff x (r)G(r, r', w) = G(r, r', w) (5) 



OO / , 



ff x (r) = -V x 



x +~ Xe(r,w), 
(6) 



i.e., G' satisfies the same differential equation as the 
one governing the electric field, with the current be- 
ing equal to G. Hence it can be written as a convo- 
lution of this current with the kernel G: G'(r,r',u/) = 
/ d 3 sG(r, s, (x>)G(s, r', uj) . Substitution of G' in this form 
in Eq. ((4]) and iteration lead to the desired Born series 
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FIG. 1: A dipole emitter in the vicinity of a rectangular 
plate. 



G(r, r', uj) = G(r, r', uj) + £ G fc (r, r', w), (7) 

fe=i 

G ^ r ' r ' w )=(n/ d3 ^) 

x G(r,si,w)G(si,s 2 ,w) • • -G(s fe ,r',w) . (8) 

This formal expansion for the Green tensor is valid for 
an arbitrary geometry, permittivity, and permeability of 
the macroscopic bodies. Obviously, one of the situations 
in which the Born series is particularly useful is when 
Xe is a perturbation to e and Xfj. 1S a perturbation to fl, 
thereby one makes only a small error cutting off higher- 
order terms. For such weakly dielectric and magnetic 
bodies, it is natural to choose 



e(r,u) = n(r,u) = 1 



(9) 



with xaO,^) =XAR(r,w) +ix\i(r,u), \xx(r,uj)\ < 1 
(A = e, fi) [cf. Eqs. ©]. This implies we restrict our- 
selves to frequencies far from a medium resonance. 

The Green tensor G corresponding to e(r, uj) = 
p,(r,L>) = 1 is the vacuum Green tensor 



G(r,r» 



a(u) 

3fc 2 " 



a = a{q) = - + — - 



q q* 



■I (al — bu (g) u)e 
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\, b = b(q) = -^ 



iq 



q q 
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(10) 
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(k= ui/c; u = r — r'; u = u/u, and q= ku). 

A description of quantities such as the emission pat- 
tern or the interatomic van der Waals forces requires 
the knowledge of the Green tensor of different positions, 
while a description of quantities such as the spontaneous 
decay rate of an excited atom or the atom-body van der 
Waals forces requires the knowledge of the Green tensor 
of equal positions. Substituting Eq. (fTO]) in Eq. ((SJ) and 
assuming that the position r lies outside the region oc- 
cupied by the macroscopic bodies, we derive for the first- 
order term in the Born expansion of the equal-position 



Green tensor 

Gi(r,r,w) 



l67r 2 / d 3 Sj ff x i(s) 



[a 2 I + 2 - 2ab)u(g) u]t 



(12) 



u 



— s; q — ku; a = a(q); b = b(q); H x i(r) = — V x 

Xp(r,w)V x +^Xs{r,u) - linear part of H x , Eq. ©]. 
Higher-order terms can easily be obtained by repeatedly 
using Eq. JTDJl in Eq. ©. 

Our system consists of an excited two-level atom sur- 
rounded by macroscopic media, which can be absorbing 
and dispersing. In the electric-dipole and rotating-wave 
approximations, the atomic decay rate reads as 



2fc| 

he 



d A Im G(r A , r A , w A )d A , 



(13) 



where d A and cj a are the atomic dipole and shifted tran- 
sition frequency, fe A = w A /c, and G is the Green tensor 
describing the surrounding media. 

In accordance with the linear Born expansion, Eqs. 
(Ill), (fl0|) . and (fl"2"l) yield, for a purely electric material, 
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u z [z 



(14) 



[Tq = fc A d A /(37rfe ) - free-space decay rate, s = (x, y, z), 
u = |s — r A |, q — fc A it, a — a(q), b = b(q)j for 
x-(z-)oriented dipole moments. Equations (jT¥J) are our 
main working equations. Just like the Born expansion of 
the Green tensor, they hold for an arbitrary geometry of 
the surrounding environment. 

Next let us be specific about the shape of the macro- 
scopic bodies. We consider a rectangular plate of dimen- 
sions d x , d y , and d z and choose a Cartesian coordinate 
system such that its origin is located at the center of 
one surface of the plate, as sketched in Fig. [TJ Then 
r" represents the spontaneous decay rate of a dipole mo- 
ment parallel to a plate surface and T ± - that of a dipole 
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FIG. 2: Atom-surface distance dependence of the normal- 
ized spontaneous decay rate of an excited atom positioned at 
(0, 0, 2a) near an infinite planar plate (solid line) and a rectan- 
gular plate (d x = d y = 10 Aa, dashed line) of equal thickness 
d z = 0.2Aa and equal Xe = XeR + *10~ 8 . Case (a) is for a 
x-oriented dipole moment, while case (b) is for a z-oriented 
dipole moment. 



moment normal to the same surface. Since no further an- 
alytical calculation in Eqs. fT4")) seems possible, we resort 
to numerical computation. 

In Fig. [2] we present the spontaneous decay rate in 
accordance with the linear Born expansion (|14p . as a 
function of the atom-surface distance for two different 
values of the permittivity. The same quantity but for 
an atom placed near an infinitely extended planar slab is 
plotted using the Green tensor given in Ref. Q. It can 
be seen that when the lateral dimensions of the rectan- 
gular plate are sufficiently large and the absolute value 
of the permittivity is sufficiently close to one (the case 
of Xe = 0.1 + il0~ 8 in the figure), the two results al- 
most coincide for both dipole moment orientations. As 
XeR increases, the agreement worsens but is still quite 
good at XsK = 0.5. Further numerical calculations indi- 
cate that in the range of |x e (o;a)| ~ 0.5, the spontaneous 
decay rate can be well approximated by the linear Born 
expansion. In Fig. [21 the atom has been moved along the 
z-axis with xa = VA = 0. When the atom is moved along 
other lines parallel to the z-axis but nearer to the border, 
the ZA-dependence of the normalized spontaneous decay 
rates behaves in a similar way as in Fig. [5] but with its 
value being generally closer to one. Having determined 
the range of permittivities where the linear Born expan- 
sion provides a good approximation to the spontaneous 
decay rate, we move on now to investigate when the plate 
can be regarded as an infinite slab. 

In Fig. [31 we gradually reduce the lateral sizes of the 
rectangular plate while keeping its thickness constant, 



FIG. 3: The same as in Fig. [2] but for different sizes of 
the rectangular plate: d x = d y — 3Aa (dashed line), 0.4Aa 
(dotted line), and 0.2Aa (dash-dotted line). In the last case, 
the plate is actually a cube. The curves for an infinite planar 
plate are shown by solid line and Xe = 0.1 + ilO -8 . 




FIG. 4: Plate-thickness dependence of the spontaneous de- 
cay rate of an excited atom located near an infinitely extended 
planar plate (solid line) and a rectangular plate (d x — d y = 
IOAa, dashed line) of equal Xs = 0.1 + il0~ 8 . The atomic po- 
sition is fixed at (0,0,0.2Aa) in the main figure, and (0,0, 5Aa) 
in the inset. 



and compare the resulting spontaneous decay rates with 
that for an infinite slab. To be on the conservative side, 
the permittivity is set equal to £(u>a) = 1-1 + ilO -8 - a 
value which is very close to one (cf. Fig. fj]). For lengths 
of the lateral sides comparable to the atomic transition 
wavelength, the infinite slab approximation starts to dif- 
fer noticeably from the linear Born expansion (see the 
figure, case of d x = d y = 3Aa), which in this situation 
is regarded as a good and nondegrading approximation. 
As the lateral sizes of the plate decrease further and be- 
come smaller than Aa, the infinite-slab approximation 
fails completely (see cases of d x = d y = 0.4Aa and 0.2Aa 
in the figure). In other words, while a rectangular plate 
with lateral sizes much larger than the atomic transi- 
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FIG. 5: Effects of the presence of a boundary in the x- 
direction on the spontaneous decay rate of an excited atom 
located near a rectangular plate of permittivity e(uja) = 
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1.5 + ilO and dimensions d z — 0.2Aa and d x 
The atom is located at (xa,0, O.OIAa). 
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method, one obtains from Eq. fT4 
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[q z = kx{z + za)]- The integrals over x and y in Eq. 1|15|) 
contain Fresnel integrals. In the limit of an infinite plate 
d x , d v — * oo, Eq. (fT5")) becomes 



^.l + ^lm 

r 16 



Xe dza 2 {q z )q z e 2 ^{l + zf 

JO 



(16) 



tion wavelength can be more or less treated as an infinite 
slab, care should be taken when the sizes are reduced 
to about or below a wavelength. This happens for both 
normal and parallel to the surface dipole moment ori- 
entations. When the rectangular plate can be roughly 
regarded as an infinite slab, the infinite-slab curve and 
the linear Born expansion curve agree better when the 
atom is placed closer to the surface (see dashed and solid 
curves in the figure). This can be explained by that the 
closer to the surface the atom is situated, the more it is 
inclined to see the plate as infinite. 

Besides the dependence on the lateral sizes, whether 
a rectangular plate can be treated as an infinitely ex- 
tended one clearly depends on its thickness as well. It 
is intuitively obvious that even when plate lateral sizes 
are much larger than the atomic transition wavelength, 
the plate cannot be treated as extending to infinity if its 
thickness is comparable with the lateral sizes. In Fig. 2J 
we compare the c^-dependence of the spontaneous de- 
cay rate for a rectangular plate with that for an infinite 
slab. The agreement between the two curves, being very 
good for sufficiently thin plates, gradually worsens with 
an increasing plate thickness. The disagreement is al- 
ready noticeable at d z ~ Aa - a value which is still much 
smaller than the lateral sizes d x — d y — IOAa, and it 
sets in earlier for Y 1 - than for Y". The two calculations 
predict quite different large-thickness limits. This means 
that in the case of thick plates, one must take into ac- 
count the boundary conditions in the lateral directions 
in order to obtain reliable results. 

Account of the boundary conditions in the lateral di- 
rections also gives rise to some curious beating in the 
d z -dependence of the spontaneous decay rate, which is 
especially visible when the atom has a dipole moment 
oriented parallel to the surface and is situated somewhat 
away from the surface (see Fig. |H inset). Let's have 
a closer look at, say, the decay rate for a dipole ori- 
ented parallel to the surface. Using the stationary phase 



As a function of z/Xa, the integrand in Eq. p^|) has a 
period of -z. Since this period is ^-independent, the re- 
sulting integral must exhibit oscillations with the same 
period, as confirmed by Fig. 0] solid curves. These oscil- 
lations survive for plates of finite lateral extensions (Fig. 
HI dashed curves). The beating clearly arises from the 
finite values of d x and d y . Note that as a function of z, 
the inner integrands in Eq. (|15p have a 'period' that is 
z- dependent and that increases with increasing z. 

Next we turn to elucidating the influence of the bound- 
aries in the x- and y-directions on the spontaneous decay 
rates. As can be seen from Fig. [5j where an edge is 
present at xa = 5Aa, the decay rates exhibit oscillations 
near the boundary with a particularly strong magnitude 
right on either side of it, and damping tails. The oscilla- 
tions are more pronounced for a dipole moment oriented 
parallel to the (xj/)-plane than for a z-oriented dipole 
moment. One can notice that when the projection of the 
atomic position on the rry-plane lies outside and suffi- 
ciently far from the boundaries, the spontaneous decay 
rate approaches one in free space, as it should. 

In summary, using the Born expansion of the Green 
tensor, we have considered the decay rate of an atom lo- 
cated near a plate of rectangular shape. We have shown 
that a rectangular plate can be treated as extending to in- 
finity only when its lateral sizes are much larger than the 
atomic transition wavelength, its thickness sufficiently 
small, and the atom is located close enough to the plate 
surface. We have also shown that a boundary in the lat- 
eral directions can give rise to significant modifications 
of the decay rate in either side of it. The first-order Born 
expansion remains quite reliable even at a value of per- 
mittivity as high as 1.5. Inclusion of higher-order terms 
would allow one to investigate more dense media. 
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